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ABSTRACT 

Many of the currently available equations of state for core-collapse supernova simulations give 
large neutron star radii and do not provide large enough neutron star masses, both of which are 
inconsistent with some recent neutron star observations. In addition, one of the critical uncertainties 
in the nucleon-nucleon interaction, the nuclear symmetry energy, is not fully explored by the currently 
available equations of state. In this article, we construct two new equations of state which match 
recent neutron star observations and provide more flexibility in studying the dependence on nuclear 
matter properties. The equations of state are also provided in tabular form, covering a wide range in 
density, temperature and asymmetry, suitable for astrophysical simulations. These new equations of 
state are implemented into our spherically symmetric core-collapse supernova model, which is based 
on general rclativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. The 
results are compared with commonly used equations of state in supernova simulations of 15 and 40 M Q 
progenitors. We do not find any simple correlations between individual nuclear matter properties at 
saturation and the outcome of these simulations. However, the new equations of state lead to the 
most compact neutron stars among the relativistic mean-field models which we considered. The new 
models also obey the previously observed correlation between the time to black hole formation and 
the maximum mass of an s = 4 neutron star. 



1. INTRODUCTION 

Core-collapse supernovae (SNe) are some of the most 
energetic events in the Universe. They release several 
times 10 53 erg in neutrinos, the gravitational binding 
energy difference between iron core and neutron star, 
and in case of SN explosions several 10 51 erg kinetic en- 
ergy of the ejected material. The latter are related to 
the standing accretion shock revival, which forms when 
the super-sonically collapsing stellar core reaches nuclear 
matter densities and bounces back. The shock wave ini- 
tially propagates out of the stellar core and thereby loses 
energy due to the dissociation of infalling heavy nuclei 
and electron neutrino escapes when it crosses the neu- 
trinospheres. Several explosion mec hanisms have been 
discus sed, i.e magneto-rotational by iLeBlanc k Wilson 
197C ) , dissipation of sound waves bv iBurrows et a!* 
2006f ). an d the standard s cenar io driven by neutrino 



heating bv lBethe k Wilson! JT9p). 

The equation of state (EOS) is one of the critical and 
highly uncertain microphysical inputs for modeling core- 
collapse supernovae. The EOS in SN simulations has to 
handle several intrinsically different regimes. For tem- 
peratures below about 0.5 MeV, time-dependent nuclear 
reactions are important in determining the nuclear com- 
position. Above 0.5 MeV, nuclear statistical equilibrium 
can be applied, where the dependence of the EOS can be 
reduced to the three independent variables: temperature 
T, baryon number density Ub and proton-to-baryon ra- 



tio (or equivalently the electron fraction) Y e . Finally, at 
densities close to and above nuclear matter density, the 
transition to a state of matter composed of deconfined 
quarks may take place. 

While accurate neutron star mass measurements of 
double pulsar systems are plentiful, reliable radius 
measurement s have only recently become available. 
I Steiner et al.l (|2010ft obtained quantitative constraints 
on the EOS of nuclear matter from mass and radius 
measurements from quiescent low-mass X-ray binaries 
and objects with photosphcric radius expansion bursts. 
Current neutron star radius observations suggest small 
radii and lower pressures just above the nuclear satura- 
tion density. Small neutron star radii could be interest- 
ing for the supernova dynamics, because expects more 
gravitational binding energy. Unfortunately very few 
available supernova EOS give (i) nuclear matter prop- 
erties consistent with those inferred from experiment, 

(ii) maximum neutron star masses large enough to be 
consistent with the recent measurement of a 1. 97 so- 
lar mass neutron star in iDemorest et al.l (|2010l ) . and 

(iii) masses and radii c onsistent wit h rece nt neutron 
star observations of e.g. iSteiner et al.l (|2010f ). Further- 
more, most available EOS tables are based on inter- 
actions with larger values of the density derivative of 
the nuclear symmetry energy, L, even though lower val- 
ues are su ggested by both int ermediate-energy heavy-ion 
collisions (Tsang et al. 2009), chiral effective field the- 



ory (iHebeler et al.ll2010t iTews et all 120121) and neutro n 
star radii (|Steiner fe Gandolfill2012t ISteiner et al.ll2012l) . 

We should note that there are potential systematic un- 
certainties in the neutron star mass and radius measure- 
ments which are not yet taken into account. The re- 
lationship between the Eddington flux and the point at 
which the photosphe re returns to the ne utron star surface 
is not under control (jSteiner et al.ll2010l ). the evolution of 
the spectr um during the tail of th e burst is not well un- 
derstood (|Suleimanov et aH 120111 ). and the value of the 
factor which corrects for the fact that the X-ray spec- 
trum is not a black body may also modify inferred radii. 
Also, the X-ray spectrum may be modified by accretion 
and violati ons of the assumed spherical symmetry. Re- 
cent work (|Steiner et al.l I2012D finds that neutron star 
radii may be as large as 13 km. 

In this article, we construct new EOSs which are 
both consistent with experimental nuclear data and 
recent neutron star mass observations. We parame- 
terize nucleonic matter with a new relativistic mean 
field (RMF) model and describe nuclei and non- 
uniform nuclear matter with the s tatistical model from 
iHempel fe Schaffner-Biehchl polof l. The final EOS is 
also provided in tabular form, covering a wide range 
in density, temperature, and electron fractions. We ap- 
ply these EOS tables in core-collapse SN simulations of 
intermediate- and high-mass progenitors. 

2. EQUATIONS OF STATE 

The most com m only used EOSs are that of 
lLattimer fe Swestvl (|1991|) (hereafter LS), which is 
based on the compressible l iquid-drop model includ- 
ing surface effects, and H. iShen et all (|1998|) (here- 
after S TOS), which is based o n the TM1 RMF inter- 
action (jSugahara fe Tokil Il99l and uses the Thomas- 
Fermi approximation to describe non-uniform nuclear 
matter. Both EOSs simplify the baryon composi- 
tion using the single-nucleus approximation (SNA) for 
heavy nuclei, and ignore all light nuclei except for 
alpha particles. There are several studies of the 
differences between thes e two EOSs in core-c o llapse 
simulations (see, e.g.. ISumivoshi et al.l [20061 120071: 
Fischer et ail^OOa lO^Connor fc Ottll20TlI IHempel etaD 



TABLE 1 
RMF Model parameters 



2012lh " More recently, several new EOS based on 
RMF interactions have become available. The ne w 
hadronic EOS table s of G. Shen (IShen et al.l [2^TTblial) 
are based on NL3 (lLalazissis et al l 11997ft and FSUg- 
old (|Todd-Rutel fc Piekarewiczl 120051 ) RMF interactions 
with nuclei des cribed in the Hartree appoximation . The 
EOS model o f IHempel fc Schaffner-Bielichl pOlOh and 
IHempel et all ([20120 (hereafter HS) is based on the sta- 
tistical approach and is also used in the present study. 
It goes beyond SNA by including the detailed nuclear 
composition, based on experimentally measured nuclear 
masses as well as different theoretical mass models. Ta- 
bles are available for seven different RMF pa rameteriza- 
tions, including TM1, FSUGold, and TMA (jToki et al.l 
[1995] ) and the new parameterizations SFHo and SFHx 
developed in the present article. 

In addition to the hadronic EOS a first-order 
phase transi t ion to qua rk matter has been studied in 
iSagert et all (|2009D and iFischer et all (|2011l ). where it 
has been demonstrated that it can trigger explosions even 
in spherically symmetric supernova models. Note that 
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the hybrid EOSs used result in extremely compact neu- 
tron stars. All of these EOS studies in SN simulations 
explore also very massive progenitor stars which in the 
end collapse to a black hole, for which the post-bounce 
time until that moment can be used as an observable to 
constrain characteristics of the phase transition. 

There are several critical parameters for character- 
izing the equation of state of hadronic matter, and 
some of the most relevant and yet uncertain parame- 
ters are the nuclear incompressibility, K, the symmetry 
energy at the saturation density, J, and the logarith- 
mic derivative of the symmetry energy L. The com- 
pressibility was an i mportant paramete r for early core- 
collapse simulations (jBaron et al.l 119851 ). and the com- 
pressibilities of the LS EOS tables are available for 180, 
200, and 375 McV. However, the lowest and highest 
of these values are far outside the curre ntly acceptable 
range of 240±10 MeV (jColo et al.l 120041) . We note that 
there is still some model d ependence in extract i ng thi s 
value iPiekarewiczl (|2010l ): iKhan fc Margueronl (|2012f ). 
Early simulations also suggested that the symmetry en- 
ergy is important (jSumivoshi et al.lll995l ). Recent con- 
straints on J suggest 28 < S < 34 from a combina- 
tion of constraints from experiments, theo ry, and obser- 
vations of neutron star masses a nd radii dHebeler et al.l 
20101 ISteiner fc Gandolfil [20ll lLatt mer fc Liml l2012t 

t ai.il: 



Tsang et al J l2012t ISteiner et alJ[20rl ITews et al^oTl ".' 

These works also imply constraints on L, but these con- 
straints are more model-dependent and are not always 
consistent with each other. Only one of the original su- 
pernova EOS tables, the LS table with K = 220 MeV, 
obeys the current constraints on these EOS parameters. 

For the new EOS SFHo and SFHx, we choose to use a 
covariant Lagrangian based on the Walecka model where 
nucleons interact via the exchange of a, w, and p mesons 
in the mean-field approximation. The non-linear Walecka 
model is well-known to have a limited range of variation 
in the isospin sector and the nuclear symmetry energy 
is controlled entirely through the coupling of the nucle- 
ons to the p meson. In order to provide more flexibil- 
ity, several additional terms like p 4 and a 2 p 2 have been 
consid ered. We use the parameterization in lSteiner et al.l 
(|2005|) which provides enough freedom to modify the low- 
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TABLE 2 

Properties at saturation density and neutron star properties for the the different EOSs under investigation. The 

definition of all the quantities is given in the text. 
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TABLE 3 

Properties of nuclei for SFHo and SFHx 





E (208p b) 
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[MeV] 
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SFHo 


-7.76 


5.44 


-8.60 


4.19 


SFHx 


-7.87 


5.41 


-8.55 


4.19 
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and high-density parts of the isospin sector separately. 
There are 17 parameters including: the scalar-isoscalar 
meson mass m CT , the standard non- linear Walecka model 
couplings g a , g u , g p , b, and c, the two fourth-order vector 
meson couplings £ and £, and the parameters ai_6 and 
6i_3 which control the symmetry energy as a function of 
density. These parameters are varied to ensure that both 
EOSs have saturation properties which agree with that 
predicted from nuclear masses and giant monopole res- 
onances. In addition, both predict binding energies and 
charge radii for 208 Pb and 90 Zr that are within 2 % of 
the experimental values. We ensure that the pressure of 
neutron matter is always positive and always increases as 
a function of the density. We ensure that the maximum 
mass is la rger than 1.93 solar m asses (the lower 1 — a 
limit from lDemorest et all ()2010t) ). We compare our re- 
sults to LS(180) and FSUGold even though these are 
inconsistent with this maximum mass because they are 
still commonly used in the supernova community. The 
requirement that the speed of sound is not superluminal 
is automatically enforced in this fully covariant model. 
In our baseline model, SFH o, we also fit the mo st prob- 
able mass-radius curve from lSteiner et all (|201Q[ ). and in 
our extreme model, SFHx, we attempt to minimize the 
radius of low- mass neutron stars yet remaining consistent 
with the other constraints given above. This forces the 
value of L for SFHx to be on the lower edge of the typi- 
cal range of 20-120 MeV. We always take the mass of the 
neutron to be 939.565346 MeV, the mass of the proton to 
be 938.272013 MeV, the mass of the w meson to be 762.5 
MeV, the mass of the p meson to be 770 MeV, a nd we 
use he = 197.3269631 MeV fm (|Mohr et all 12008ft . The 
full parameter list for both models is given in Table [TJ 
The saturation properties, nuclear binding energies and 
nuclear charge radii for our models are given in Tables [2] 
and [3] We remark that the values given in Table [2] can 
be slightly different to previously published ones, e.g. by 
iHempel et al.l ()2011[ ). due to a different treatment of the 
nucleon rest masses. In the present work we are using 
the measured masses from above for all of the HS EOS, 
to obtain the correct low-density limit. This treatment 
also leads to a slight splitting of the neutron and proton 
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Fig. 1. — Pressure of zero temperature beta-equilibrium matter 
as a function of the total energy density. The gray region give s 
the 1 <r confidence limits on the EOS from ISteiner et al.l (2010). 
The curve labeled TM1 applies for both the HS(TMl) and STOS 
EOSs. 

effective masses, m* and m* respectively. The value of 
J for LS in Table p differs from the published value by 
lLattimer fc SwestvT(| 19911 ) of 29.3 MeV. These authors 
computed J as the energy difference between neutron 
and nuclear matter whereas we are calculating J as the 
second derivative with respect to Y e at the saturation 
point. In Table [21 we also give the saturation density n° B , 
the binding energy of nuclear matter Eq , the skewncss of 
nuclear matter K' , the reduced neutron and proton ef- 
fective masses at saturation, the radii of 1.4 solar mass 
neutron stars at T = 0, the maximum mass at T = 0, and 
the maximum mass for stars with constant entropy s = 4 
in beta-equilibrium without neutrinos. The T = EOSs 
for all of the models constructed with nucleon degrees of 
freed om are given i n Figu red! along with the constraints 
from lSteiner et al.l (120101 ). Figure^ gives the correspond- 
ing mass-versus radius c urves and also the constraints 

1 ipolnh . 



from lSteiner et al.1 

For densities belo w saturation density we apply t he sta - 
tistical model from IHempel fc S chaffnc r-Bielichl (|2010l ) 
(HS) for the description of non- uniform nuclear matter in 
nuclear statistical equilibrium (NSE), i.e. for the forma- 



4 




Fig. 2. — Neutron star mass-radius curve for modern supernova 
equations of state. The red (green) region o utlines the one (two) 
a confidence limits from[Stcincr et al. (201(3). 



tion of light and heavy nuclei within the gas of unbound 
nucleons. For the unbound nucleons, we utilize the SFHo 
and SFHx relativistic mean-field interactions. At low 
densities, the description of nu clei is based on m easured 
experimental binding energies (|Audi et al.l 12001 . which 
are combined with theoretical nuclear structure calcu- 
lations for exotic nuclei withou t measured data. Here 
the finite range droplet model of lMoller et~al"1 (|1995l ) was 
chosen because of its excellent reproduction of experi- 
mental binding energies, with a rms deviation of only 
0.669 McV. Due to the use of nuclear structure data, 
shell effects are automatically included. HS goes beyond 
the single nucleus approximation and utilizes a distribu- 
tion of different nuclear species, and the results for light 
nuclei are in agreement with more sophistic ated quan- 
tum many-body models (iHempel et al.ll201ll). Also, the 
recent experimental study of iQin et all (|2012[ ) indicates 
that the HS model is well suited for the description of 
matter at finite temperature and densities around a few 
tenths of saturation density. At even larger densities in 
the HS model the disappearance of nuclei and smooth 
transition to uniform nuclear matter is assured by an 
excluded volume description. Finally wc calculate the 
EOS in tabular form, covering densities from 10~ 12 to 10 
fm~ 3 , temperatures from 0.1 to 160 MeV, and electron 
fractions from to 0.6, including detailed information 
about the nuclear composition and the thermodynamic 
properties. The tables are suitable for use in astrophysi- 
cal simulations and are available online 

The resulting EOS will be compared below with the 
LS EOS, with the different compressibilities 180 MeV 
(LS180), 220 MeV (LS220), and 375 MeV (LS375), and 
with STOS. Moreover, we will also include into our com- 
parison r esults obtained using the quark-hadron hybrid 
EOS from Fis cher et al.1 (|2011| ). We select the model with 



bag constant B 1 ' 4 = 155 MeV and including corrections 
from the strong coupling constant, as = 0.3 (hereafter 
QB155a!503), where the phase transition to quark mat- 
ter takes place at nuclear saturation density for temper- 
atures around 10 MeV and Y e ~ 0.3 (for details, see 
IFischer et all 1201 If ). The hadronic part of this EOS ta- 
ble is based on the STOS EOS, and these two EOS are 
identical at sub-saturation densities where quarks are not 
present. 

3. CORE-COLLAPE SUPERNOVA SIMULATIONS 

In this section, we will compare results from SN sim- 
ulations obtained using the SFHo EOS with the stan- 
dard EOS LS180 and the two TM1 RMF parameteriza- 
tions STOS and HS. Furthermore, we will also compare 
SFHo with the hybrid EOS QB155a s 03, for which ex- 
plosions were obtai ned recently even in spherically sym- 
metric simulations (jFischer et aLl[20TTl ) . 



TABLE 4 

Neutrino reactions considered including references. 



Rcaction a 



References 



v e + n -> p + e 
P e + P — ► n + e+ 
(A, Z - 1) -> (A, Z) - 



v + N 
(A,Z) 

i/ -X- 



.JV- 



u' + N 
v> + {A,Z) 



Bruenn (1985) 
Bruenn (1985) 
Jjaiiganke^et^al^ (2003), 
Hix et al. (2003) 
Bruenn (1985) 
Bruenn (1985) 
Bruenn f 1985) . 



Mczzacappa 
Mczzacappa 



-N - 



N + N 



^ezzacamoa 



Bruenn (1985) 
Messer 



Hannestad fc Raffelt (1998) 
Buras et al. (2003) 



'Note: v = {u e , u e , v^/t, u ^/t} an d N = { n iP} 



3.1. Supernova model 

Our corc-collapsc SN model, AGILE-BOLTZTRAN, 
is based on general relativistic radiation hydrodynam- 
ics in spherical symmetry. It e mploys three-flavor Boltz- 
mann neutrino transport (see iLiebendoerfer et al.l 120041 
and references t herein) . We us e the standard weak pro- 
cesses following iBruenr] (|1985f ). see Table |4] for details. 
In addition, we include the impr oved rates for el e ctron - 
capt ures on heavy n uclei from iLanganke et al.l ()2003l ) 
and iHix et akl (I2003D , weak m agnetism and nucleon re- 



1 See http : //phys-merger .physik .unibas . ch/~hempel/eos .html. 



coil based on lHorowit3 (I2002D . and the annihilation of 
trapped electron neu trino pair s has been impleme nted in 
IFischer et all (|2009D following iBuras et al.l (|2003[) . 

For NSE conditions (T > 0.45 MeV), we implement 
the baryon EOS tables specified above. For non-NSE, 
we assume the ideal gas of 28 Si for the baryon EOS. On 
top of the baryons, also for NSE, contributions from elec- 
trons, positrons and phot ons are added to the EOS using 
iTimmes fc Arnettl (|1999l ). Recently, this Si-gas approxi- 
mation has been replaced by a nuclear reaction network, 
based on the nuclear composition given by the progeni- 
tor model. It allows, e.g., for a smooth NSE-to-non-NSE 
transition as well as to simula te a large domain o f the 
progenitor star (for details, see IFischer et al.ll2010D . 

The simulations wc will discuss further below arc 
launched from iron-core progenitors. We use the 
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15 M Q model from IWooslev et all ([2002) for regu- 
lar core-collapse supern ovae and the 40 M Q from 
IWooslev fc Weaver] (|1995l ) for the black-hole formation 
scenario. None of the spherically symmetric simulations 
using the purely hadronic EOS results in an explosion for 
the considered simulation times. 

3.2. Simulation results 

A detailed supernova EOS comparison study for LS180 
and STOS as well as several HS EOS tables with different 
R MF parameteri z ation s has been published very recently 
in iHempel et al.l (|201 2T ) . Here, we extend their analysis 
and include in addition SFHo/SFHx. 

During the core-collapse phase, the composition is 
dominated by heavy nuclei. At low temperatures in non- 
NSE the EOS based on the ideal gas of Si-nuclei has been 
applied in all simulations except for HS(TMl) where we 
used our simplified nuclear reaction network. Hence, all 
EOS lead to very similar structures in the simulations at 
low temperatures, see e.g. the equal entropy per baryon 
profiles in the outer layers (Mb > 1.5 M Q ) as illustrated 
in Fig. [3] (except HS(TMl)). Note also that LS(180) has 
a different non-NSE treatment based on the ideal Si-gas. 
Moreover, structure differences arise at higher tempera- 
tures in NSE, which will be discussed i+n the following 
paragraph. In general, low-density differences (i.e. on 
top of the bounce shock 0.7 < M B < 1.5 M in Fig. [3j 
between STOS and HS(TMl) are related to a different 
description of heavy nuclei. The statistical approach 
of HS is largely based on experimentally known masses 
including shell effects. It resembles the ideal gas of 
56 Ni/ 56 Fe at the transition to non-NSE by construction. 
In contrast, the Thomas-Fermi approximation of STOS 
gives heavy nuclei which are too strong ly bound and 
does n ot perform well at this transition. IHempel et al.l 
(|2012[) explained how these differences affect the Y e evo- 
lution. Furthermore, entropy differences between STOS 
and HS(TMl) originate from the missing kinetic entropy 
contribution of heavy nuclei in STOS. Since we use the 
same low-density description for nuclei as in HS, simu- 
lations using SFHo and SFHx result in the same con- 
ditions during the core-collapse phase. The simulation 
results for QB155as03 are similar to STOS at low den- 
sities where quarks are not present. Even though they 
are based on the same low-density EOS, small differences 
have emerged because the core has already entered the 
quark-hadron mixed phase. 

Only a few tens of milliseconds before core bounce, 
nuclear saturation density is reached at the center. At 
these conditions high-density EOS differences become 
large which in turn lead to different dynamical evolu- 
tions. We define the moment of core bounce when the 
central density reaches its maximum before shock break- 
out. The EOS with highest central density, LS(180), 
results in the smallest enclosed baryon mass (see Fig. [3]). 
The two TM1 EOSs, STOS and HS(TMl), result in the 
lowest central densities. Differences between STOS and 
HS(TMl) are due to the different description of heavy 
nuclei, as already discussed above. Furthermore, the 
presence of light nuclei, which are taken into account 
in HS, also has an impact on thermodynamic properties. 
Comparing STOS and QB155cts03, the most pronounced 
differences arise above nuclear saturation densities which 
are due to the presence of deconfined quarks (in partic- 




FlG. 3. — Bounce profiles of selected quantities for the 15 Mq 
models, comparing the different EOS under investigation. 



ular s-quarks). They soften the high-density EOS and 
result in higher central densities and signi ficantly lower 
centr al Y e at core bounce (for details, see iFischer et al.l 

HqIH). 

One of our new models, SFHx, gives the lowest elec- 
tron fractions for M < 0.25 Mq, due to its low slope of 
the symmetry energy L. Compared with LS(180), SFHo 
and SFHx reach similarly high central densities and low 
central Y e at core bounce (see Fig. [3] and further discus- 
sion below). Note that LS(180) has lowest incomprcss- 
ibility K and symmetry energy J among all EOS under 
consideration. However, the baryon mass enclosed inside 
the bounce shock for the two new EOS differs by about 
0.1 M and is actually more similar to the TM1 results, 
which has an extremely high K and J. These results 
provide no clear correlations between individual nuclear 
matter properties of the EOS (see Tabled and the dif- 
ferent conditions, e.g., central density and Y e obtained 
at core bounce for these EOS (sec Fig. UJ. It seems dif- 
ficult to disentangle individual nuclear EOS properties 
from the different conditions obtained in simulations of 
stellar core collapse. 

The central object formed at core bounce, i.e. the pro- 
toneutron star (PNS), is hot and lepton rich in which 
sense it differs from the final remnant neutron star. After 
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bounce, the shock starts to propagate outwards with ini- 
tially positive velocities. Simultaneously, mass accretion 
from the outer part of the stellar core continuously grows 
the mass of the PNS. Moreover, the expanding bounce 
shock loses energy due to heavy-nuclei dissociation and 
neutrino emission. The neutrino emission is related to a 
large number of electron captures on protons during the 
shock passage across the neutrinospheres. It releases a 
burst of v e of several 10 53 erg s _1 for a short timescale 
between 5-20 ms after core bounce. Both sources of 
energy loss turn the expanding and dynamic (i.e. ac- 
companied with matter outflow) bounce shock into an 
accretion front, the standing accretion shock (SAS), al- 
ready between 5-10 ms after core bounce. The later 
PNS evolution is determined from mass accretion and 
the subsequent PNS compression which leads to continu- 
ously rising central density and temperature. For a given 
mass accretion rate, determined from the progenitor, and 
otherwise identical simulation setup, the timescale for 
the PNS compression is directly related to the EOS. On 
timescales on the order of several 100 ms, EOS differ- 
ences lead to different neutrinospheres and shock radii, 
displayed in Fig. 01 The neutrinospheres can be used to 
characterize the size of the central protoneutron star as 
they are located in its outer envelope. 

At 100 ms post bounce, LS(180) has the smallest shock 
radius at 130 km because the shock had lower energies at 
bounce. However, its contraction proceeds initially, up to 
about 200 ms post bounce, somewhat slower than those 
for the EOS HS(TMl) and SFHo/SFHx. This behavior is 
related to the rapid PNS contraction obtained for the lat- 
ter. By comparing with STOS, one can identify the effect 
of the different treatment of heavy nuclei and in addition 
the consistent inclusion of lig ht nuclei for HS and also for 
SFHo/SFHx (for details, see lHempel et al.H2012t ). These 
differences of the low-density EOS result in an initially 
accelerated PNS contraction of HS(TMl). STOS leads to 
the generally slowest PNS contraction and hence slowest 
shock contraction as well. Differences between STOS and 
QB155as03 arc related to the presence of quarks, in par- 
ticular strange quarks, at and above saturation densities. 
These become abundant during the PNS contraction as 
the PNS domain occupied by quark matter grows con- 
tinuously during the accretion phas e (for details about 
the scenario, see iFischer et al.ll2011|) . The high-density 
EOS is dominated by the quark matter description after 
about 200 ms post bounce. This results in an accelerated 
contraction of the PNS for QB155as03 in comparison to 
STOS (see Fig. [J). Note that at about 300 ms post 
bounce, the PNS becomes gravitationally unstable and 
collapses. A massive quark core forms at the PNS cen- 
ter, which causes the collapse to halt and an additional 
shock wave forms at saturation density. The subsequent 
shock expansion and even acceleration, accompanied by 
mass outflow, determines the onset of explosion for this 
model. Since none of the purely hadronic EOSs lead to 
an explosion, we only show and discuss results up to the 
moment of PNS collapse for QB155as03. 

It is difficult to predict astrophysical simulation results 
from only nuclear matter properties at saturation density 
and zero temperature. For example both HS(TMl) and 
STOS arc based on the same nucleon interactions TM1 
(see Table HI). However, the additional inclusion of light 
nuclei and the different description of heavy nuclei, rel- 
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evant for the low- and intermediate-density EOS have a 
significant impact on the conditions at bounce and then 
also on the later PNS structure and contraction behavior. 
The early PNS evolution for our new parametcrizations 
SFHo and SFHx, up to about 50 ms post bounce, re- 
sembles those of HS(TMl) (see the neutrinospheres in 
Fig. HJ). Later, the PNS contraction proceeds even faster 
than for HS(TMl), reflecting the softer high-density EOS 
of SFHo/SFHx. On timescales of several 100 ms, the 
PNS and shock contractions slow down for all the RMF 
EOSs as a consequence of the continuously growing en- 
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closed mass. This effect is significantly weaker for the 
non-rclativistic LS(180) EOS. 

Fig.[5]shows the evolution of the compactness, the ratio 
of the enclosed gravitational mass to the radius. We eval- 
uate these quantities at the electron ncutrinosphere, rep- 
resentative of the protoneutron star surface. The EOSs 
which result in the most compact protoneutron stars at 
core bounce and few tens of milliseconds after are those 
of SFHo and SFHx, closely followed by HS(TMl). At 
core bounce and up to 50 ms post bounce, LS180 results 
in the least compact protoneutron star. However, dur- 
ing the later post-bounce evolution for LS180 the com- 
pactness of the central protoneutron star grows fastest 
among all of the EOSs under exploration. For HS(TMl) 
the protoneutron star contraction proceeds on the same 
timcscale as STOS, only on a slightly more compact 
magnitude. This can be attributed to the more com- 
pact, configuration at core bounce for HS(TMl). STOS 
and QB155as03 proceed along identical lines until about 
100 ms post bounce, when quark matter begins to domi- 
nate the protoneutron star evolution. Note, at that mo- 
ment a quark-hadron phase transition to strange quark 
matter has not been fully achieved, i.e. most mass inside 
the protoneutron star is only in the quark-hadron mixed 
phase. The EOSs with most compact protoneutron stars 
at bounce and during the early post-bounce evolution 
are SFHo and SFHx which remain the most compact 
ones until about 300 ms post bounce, when only LS180 
becomes more compact. Up to this moment, SFHo and 
SFHx have the same compactness, and differences ob- 
served during the later evolution remain small. 

In addition to the simulations using the 15 Mq progen- 
itor, we also include the mor e massive 40 M© progeni- 
tor from I Wooslev et all <|2002f ) into our EOS comparison 
in supernova simulations. This progenitor has been dis- 
cussed before within the scenario of black- hole formation, 
also explored for the LS and STOS supernova EOSs (see , 
e.g.. iSumivoshi et all [20061 [20071: iNakazato et all 120081; 



Fischer et al. l l2009HO'Connor fc Ottll2011l:lHempel etaj] 
20121) . The moment of black hole formation is defined, 
wit hin our general relativist ic framework (for details, 
see lLiebendoerfer~e t al. 2004]), as the central lapse func- 
tion approaches zero and stable numerical solutions for 
the evolution equations cannot be obtained anymore. 
The lapse function determines, e.g., time dilatations 
and gravitational Doppler shift effects. Note that core- 
collapse events which result in a black hole instead of an 
explosion are as bright as normal supernova explosions 
with respect to neutrino signal. The possible future de- 
tection of such a neutrino signal from a Galactic event 
ma y be used to fur t her co nstrain the high-density EOS, 



ma y be used to tur tner co 
see lNakazato et all (|2010h 



The moment of black-hole formation is approached 
during the post-bounce evolution, during which the mass 
enclosed inside the PNS grows continuously in the ab- 
sence of an explosion. The timescale for the central den- 
sity (and also temperature) rise is given by the PNS con- 
traction behavior, which in turn depends on the nuclear 
EOS. We list selected properties at the onset of PNS col- 
lapse in Table[5]for all hadronic EOS under investigation. 

The PNS compression can be illustrated best via the 
central density evolution, shown in Fig. [HI comparing the 
different hadronic EOSs under investigation. The values 
of the central density cluster around only two different 
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Fig. 6. — Central density evolution for the 40 Mq model, com- 
paring the hadronic EOS under investigation. 



values, 3.2xl0 14 and 2.5xl0 14 g/cm 3 , immediately after 
the bounce. It is likely that this is related to the prop- 
erties of the EOS at low-density, and thus also to the 
radii of very low-mass neutron stars. A similar (but not 
exact) clustering is observed in the lower-right hand cor- 
ner of Fig. H where STOS, HS(TMl), and HS(FSUGold) 
have larger radii for 0.4 solar mass neutron stars. This 
grouping lasts only for 300 ms after which time the EOS 
at higher densities becomes important. 

HS(TMA) has a higher incompressibility but lower 
symmetry energy than STOS (which uses TM1) and 
leads to a significantly shorter (about 300 ms) accretion 
time until black- hole formation (see Fig. [6]) . This would 
indicate that the symmetry energy plays the dominant 
role. However, SFHo and SFHx do not obey this trend. 
Even though SFHo has a higher incompressibility and 
higher symmetry energy than SFHx, it leads to a shorter 
(about 50 ms) accretion time until black-hole formation. 
Instead of using nuclear matter properties, one could ex- 
pect that our results could be explained by the maximum 
mass of cold neutron stars, which is an integrated quan- 
tity of the EOS. However, by comparing the times until 
black hole formation from Table [5]) with the numbers in 
Table [21 it turns out that there is no monotonic correla- 
tion. For example the collapse of FSUgold occurs about 
160 ms later than for LS(180), even though the maximum 
mass of LS(180) is 0.1 M larger. 

Note that the two RMF EOS parameterizations, TM1 
and TMA, fail to fulfill both zero temperature maxi- 
mum mass and radius constraints, while LS(220) is on 
the edge of acceptance, see Fig. [2j The RMF EOS based 
on TM1 and TMA result in sufficiently high maximum 
masses but give at the same time radii which are too 
large. With the black-hole formation scenario explored 
here, these large radii are due to slower PNS contraction 
which results in an extended post-bounce mass accretion 
period until black- hole formation (see Fig. [6]). On the 
other hand, FSUgold results in reasonably small neutron 
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TABLE 5 

Selected quantities at the onset of PNS collapse and 

BLACK-HOLE FORMATION. 



EOS 


[s] [10 15 g/cm 3 


[MeV] 


M B ± 

[M Q ] 


M G 2 

[M ] 


LS(180) 


0.415 


1.292 


29.978 


2.227 


2.133 


LS(220) 


0.521 


1.324 


31.446 


2.350 


2.233 


SFHo 


0.551 


1.067 


46.691 


2.477 


2.347 


HS(FSUgold) 


0.571 


1.058 


48.104 


2.465 


2.341 


SFHx 


0.625 


0.803 


40.830 


2.552 


2.424 


HS(TMA) 


0.737 


0.943 


46.708 


2.626 


2.466 


STOS 


1.028 


0.769 


49.705 


2.864 


2.652 



Notes: 

a time post bounce 
b central baryon density 
c central temperature 
d enclosed baryon mass 
c cncloscd gravitational mass 



star radii but fails to fulfill the current maximum zero 
temperature mass constrain of 1.97 ± 0.04 M©. It leads 
to a significantly shorter time until black-hole formation 
than for TMA and TM1, comparable to SFHo which has 
a larger incomprcssibility but slightly lower symmetry 
energy than FSUgold. On the other hand, SFHo and 
SFHx fulfill both mass and radius constraints by con- 
struction. It is interesting that the models which ful- 
fill the mass-radius and maximum mass constraints, i.e. 
LS(220), SFHo, and SFHx, give rather similar results for 
tBH, which seems to be constrained to 500 — 650 ms. It 
appears that neutron star radius measurements constrain 
t he time to blac k hole formation for 40 M© progenitors. 

iHempel et al.l (|2012l ) demonstrated that maximum 
neutron star mass determined from the EOS in beta- 
equilibrium at s = 4 ks /baryon was strongly correlated 



with the time until black-hole formation. This time can 
be measured with currently available neutrino detectors 
given a galactic core-collapse supernova. The associated 
mass versus radius curves and their relationship with the 
time to black hole formation is displayed in Fig.[7]for the 
hadronic EOSs under investigation. The corresponding 
values are also listed in the last column of Table Note 
that the maximum masses increase for all EOS compared 
to the T = case, (compare also Figs. [7] and [2]). Pressure 
is more sensitive than the energy density to the temper- 
ature effects because of the large nucleoli mass and thus 
larger temperatures in increase the maximum mass. If we 
compare the maximum masses of the s = 4 kg /baryon 
configurations from Tabic [2] with the time until black- 
hole formation from Table [5j we find that these quanti- 
ties have the same ordering, i.e. show an (almost) strictly 
monotonic correlation. 

4. SUMMARY AND CONCLUSIONS 

We developed two new relativistic mean-field (RMF) 
interactions, based on which we constructed the two new 
supernova EOSs, SFHo and SFHx, in such a way to ful- 
fill current constraints from neutron star mass and ra- 
dius measurements. Moreover, the new EOSs are con- 
sistent with nuclear experimental constraints on matter 
near and below the saturation density. These new EOSs 
provide more variation in the set of EOS tables which 
can be used by the core-collapse supernova community 
instead of EOSs based on TMA or TM1 which are now 
ruled out by observations. 

The new EOS were implemented in core-collapse su- 
pernova simulations of massive iron-core progenitors. 
We compared the r esults with the commonly used non- 
relativistic E OS of lLattimer fc Swestvl (|1991l ) and the 
RMF EOS of iShen et al.l (|1998t ). Moreover, we include 
into our EOS comparison also the recently introd uced 
quark-hadron hybrid EOS fromlFischer et al.l (120111) and 
the RMF EOS from IHempel fc Schaffner-Bielichl (| 20101 ) 
based on the parameterizations TM1, TMA and FSUg- 
old. We compared the different EOS during the iron- 
core collapse phase, which is dominated by heavy nu- 
clei, and confirmed al ready reported differ ences between 
these EOS (see, e.g.. IHempel et all 12012] ) . We extend 
the analysis and include SFHo/SFHx. Differences be- 
come large only slightly before and after core bounce, 
when the central density exceeds normal nuclear matter 
density. The post-bounce mass accretion phase is ideal to 
study the protoneutron star contraction behavior, which 
reflects the EOS underlying nuclear matter properties for 
a given progenitor choice. We found that the two new 
EOS which give small radii for cold neutron stars also 
lead to the most compact protoneutron stars in the first 
300 ms after bounce. However, it is not easy to disen- 
tangle the relationship between individual nuclear matter 
properties, given at zero temperature near the saturation 
density, and the outcome of our supernova simulations. 

We also find the EOS classifications 'soft' or 'stiff' mis- 
leading. Implicitly what is meant by soft or stiff is that 
the EOS has a lower or higher pressure. However, core- 
collapse supernova explore a large range of densities and 
temperatures, and an EOS which has a higher pressure at 
one density and temperature may have a lower pressure 
at another density and temperature. This is particularly 
evident with the black-hole formation time as described 
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here and in iHempel et all (|2012T >. as EOSs like SFHx 
which have a low pressure at zero temperatures near the 
saturation density, have a larger time until black hole 
formation than the other EOSs because their pressure at 
s = 4 is larger. The addition of quark degrees of freedom 
only further complicates this issue. 

Moreover, we explored the possible correlation between 
the EOS and outcome of the black-hole formation sce- 
nario in the absence of a supernova expl osion. We con- 
firm the analysis of IHempel et al.1 (|20 12f ) . where only a 
correlation between finite entropy per baryon configu- 
rations of static protoneutron stars and the time until 
black-hole formation was found. Note that the neutrino 
signals of these events are as bright as ordinary super- 
nova explosions. But at the moment of black-hole for- 
mation the signal suddenly stops which makes this time 
measurable with neutrino telescopes. The possible fu- 
ture neutrino observation of such a Galactic event will 
constrain the high-density and finite-entropy EOS sig- 
nificantly, complementary to future neutron star mass 
and radius observations. 

One conclusion from this study is that it seems to be 
difficult to understand the effect of the EOS in core- 
collapse supernovae by analyzing non-exploding models. 
To tackle the impact of the EOS on the explosion mech- 
anis m, multi-dimension al simulations are necessary, see 
e.g. iMarek et al.l (120091 ). Furthermore, it remains to be 
explored how different nuclear EOS can influence long- 
term cooling to protoneutron stars after the supernova 
explosion has been launched. Therefore, the implemen- 
tation of weak processes consistent with the nuclear EOS 
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is required (see, e.g.. lReddv et al.lll998l ri999f ) . This may 
impact the neutrino cooling timescale and also the ex- 
tent of proto neutron star convectio n as studies recently 
described bv iRoberts et al.1 ()2012l ). The nuclear EOS 
may also be important for the proton-to-baryon ratio of 
the material ejected form the protoneutron star surface 
known as neutrino-driven wind relevant for the nucle- 
osynthesis of heavy elements. Moreover, a possible EOS 
impact within the cooling of protoneutron stars on the 
emitted neutrino signal may also be of relevance for neu- 
trino oscillation studies, in particular those which explore 
collective phenomena in the presence of large neutrino 
but small matter densities. 
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